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Abstract. The electronic structure of a 2D gas subjected to a tilted magnetic field, with a 
strong component parallel to the GaAs/AlGaAs interface and a weak component oriented per- 
pendicularly, is studied theoretically. It is shown that the parallel field component modifies the 
originally circular shape of a Fermi contour while the perpendicular component drive an electron 
by the Lorentz force along a Fermi line with a cyclotron frequency given by its shape. The corre- 
sponding cyclotron effective mass is calculated self-consistently for several concentrations of 2D 
carriers as a function of the in-plane magnetic field. The possibility to detect its field-induced 
deviations from the zero field value experimentally is discussed. 
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1. Introduction 



The energy level structure of quasi-two-dimensional systems in magnetic fields tilted with respect 
to the sample plane have attracted attention for many years. As shown by Maan [fij and Merlin 
0] the electron spectra can be found analytically for a special case of parabolic quantum wells. 
This simple analytically solvable model is very useful when discussing qualitative aspects of 2D 
electron physics in quantum wells, but when the semi-quantitative comparison of results with 
experimental data is the goal of investigation, more realistic models must be used and numerical 
methods of calculation employed. 

The widely accepted approach to the electronic structure of realistic quantum wells in tilted 
magnetic fields, which relies on the perturbation theory, was reviewed by Bastard Q. It can be 
outlined as follows. The Hamiltonian is decomposed into three parts, H = Hn + Hj_ + Hint- The 
first part Hn describes the motion of an electron in z— direction under the combined influence of 
the confining potential and the in-plane component of the magnetic field. If the in-plane field is 
not too strong its effect is usually included via the perturbation theory. The second part of the 
Hamiltonian H± corresponds to the standard 2D gas subjected to the perpendicular component 
of the field and can be diagonalized analytically. In most situations the third coupling term Hint 
can be completely neglected or treated as a small perturbation. Only in certain special cases, 
e.g. when crossing of Landau levels from different subbands is important, the matrix elements 
of a coupling term must be evaluated and the Hamiltonian H diagonalized numerically jlj . 

In this paper we present a slightly different approach to the electron structure of quantum 
wells in tilted magnetic fields which is appropriate when the perpendicular component of the 
applied magnetic field is weak. Assuming for a moment B = (0, B y ,0), i.e. the field exactly 
parallel to the plane of 2D electron gas, the Hamiltonian Hn of an electron confined to the x — y 
plane by potential V con f(z) can be written as 

h\\ = rtl <*> ~ eB ^ 2 + L p y + + v ^ {z) - (1) 

Its energy spectrum is formed by subbands, momentum operators p x and p y commute with Hn 
and the eigenenergies E n (k x , k y ),n = 0, 1, . . . are the functions of the quasi-continuous wavevec- 
tors k x and k y . The 2D Fermi surfaces, or more accurately the Fermi contours, are then defined 
separately for each subband by the equation Ep = E n (k x , k y ). 

The in-plane magnetic field B y does not influence the electron motion in the y— direction 
while the electrons moving in the x— direction are slowed down or accelerated by the combined 
effect of the crossed fields B y and E z = —dV con f(z)/dz, depending on the form of V con f(z). 
Consequently, it is expected that the subband separation, the 2D density of states and the 
shape of the Fermi contour will vary with increasing B y |H|. Only recently realistic self-consistent 
calculations of electron energy spectra in heterostructures subjected to parallel magnetic fields 
have been performed ||, |(|, enabling the quantitative estimate of these quantities. 

Let the weak perpendicular component B z of the magnetic field be added to the strong 
in-plane field component B y . Then the electron dynamics in the x — y plane can be described 
semiclassically, similarly as the dynamics of electrons in metals with non-spherical Fermi surfaces. 
The corresponding quasi-classical theory was originally developed by Onsager @ and Lifshitz 
||], and nowadays it is a part of standard textbooks devoted to the solid state physics (see 
e.g. H). According to this theory electrons move in the fc-space along trajectories defined by 
intersections of a Fermi surface and planes perpendicular to the applied magnetic field. In the 
2D case there is only one such line, identical with the Fermi line itself. Once the fc-space orbit is 
known, the semiclassical theory predicts that the electron real space trajectory is of an identical 
shape except of the scale factor fi/\e\B z , and rotation by n/2. The anisotropy of the Fermi 
contour due to the strong parallel field thus leads to the deviation of the electron trajectory 
from the originally circular form. 
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The modification of the the Fermi contour shape by the in-plane magnetic field should 
manifest itself through the field dependence of the cyclotron effective mass m c . This important 
characteristics of electron energy spectra is related to the Fermi area Sp surrounded by the 
Fermi contour by a semi-classical expression 

Note that the field dependence of m c is determined by the combined influence of the in-plane 
component of the magnetic field and of the shape of the confining potential V con f(z) of the 
quantum well and, therefore, the different structures with different quantum wells can be dis- 
tinguished by measuring this quantity. 

Two methods are widely used to determine the cyclotron effective mass of 2D electrons 
confined to GaAs/AlGaAs interface: the cyclotron resonance in the infrared region of optical 
spectra and the temperature damping of Shubnikov - de Haas oscillations Q]. To get insight 
into the feasibility of such experiments, we report in this paper the results of self-consistent 
calculations of the field dependence of the cyclotron effective mass in the electron layer confined 
to the interface of the standard GaAs /AlGaAs heterostructure for several concentrations of 2D 
carriers. 



2. Electronic structure in tilted magnetic fields 

As mentioned above, non-interacting electrons mobile in the x — y plane and confined by the 
potential V con f(z) in the z-direction are considered. A tilted magnetic field B = (0, By, B z ) is 
applied to the system. The corresponding one-electron Hamiltonian reads 

H = ^-(p-eA) 2 + V conf (z) (3) 
2m 

where e and m are electron charge and effective mass, respectively. We choose the vector 
potential A in the form A = {—B z y + B y z, 0,0) which can be split into a sum A = Aj_ + An 
of vector potentials A^ = (—B z y, 0, 0) and An = (B y z, 0, 0) describing the perpendicular and 
parallel components of the magnetic field, respectively. 

To bring the above Hamiltonian to a form similar to (1) we introduce new canonical momenta 

by 



= p x ~ eA± x = p x - mu z y, 

(4) 

ITy = Py — eA±y = Py 

where lo z denotes \e\B z /m. These momenta correspond to the in-plane electron motion and the 
perpendicular component of the magnetic field B z enters Hamiltonian exclusively through them. 

While the quantum mechanics is inevitable to obtain the correct picture of an electron 
motion in the z— direction, both quantum mechanical and classical descriptions are acceptable 
for the in-plane motion when the perpendicular component of the magnetic field B z is weak. 
The reason is that the weak perpendicular component yields many occupied Landau levels below 
the Fermi energy. The Landau levels near the Fermi energy are represented by the states with 
high quantum numbers and, as it is well known, they can be treated quasi-classically as well 
as quantum-mechanically. Therefore, in the first step of our approximate description the in- 
plane canonical variables and momenta (including the part Aj_ of the vector potential) will be 
considered as classical variables which commute with the Hamiltonian and enter it as c— numbers. 



3 



Introducing further the components of the wavevector, k x and k y , by tt x = fik x and ir x = hk x , 
respectively, we obtained the one-dimensional Hamiltonian identical to (1). 

Again, the energy spectrum of this Hamiltonian is formed by subbands and the eigenenergies 
E n (k x , k y ) are the continuous functions of the wavevectors k x and k y or, equivalently, of the in- 
plane canonical momenta ir x and 7r y . Since the perpendicular component of the magnetic field 
is hidden in the wave vector component k x , the Hamiltonian describes the 2D system subjected 
only to the parallel field B = (0,B y ,0). The shape of subbands is determined by both the 
confining potential and the in-plane magnetic field and can be accurately obtained by the self- 
consistent numerical calculation. 

In the second step of our approximate treatment the single-subband effective Hamiltonian 
H e ff for the in-plane electron motion is constructed from a subband energy E n (k x ,k y ) by the 
backward substitution k x — ► n x /H, k y — ► ir y /Ti. We get 

Heff = E n (ir x ,ir y ) (5) 

and this effective Hamiltonian will be used to describe the electron motion in the x — y plane 
by standard Hamilton equations x = dH e ff /dp x , y = dH e ff /dp y , p x = —dH e ff/dx, p y = 
—dH e ff/dy. The time derivatives of coordinates define the velocity components v x = x,v y = y 
and yield expressions 

IdEnldEn (R . 

Vx ~ h dk x ' Vy ~ h dk y • {b) 

The time derivatives of momenta lead to 

hk x = -\e\B z v y , fik y = \e\B z v x . (7) 

These are the classical equations of motion for an electron subjected to the Lorentz force 
e (v x B), B = (0,0, B z ). Their solutions are in the form of orbits of all energies, among 
them the orbit corresponding to the Fermi energy is the most important. A simple geometrical 
analysis allows to derive the following properties of orbits. 

As already mentioned, an electron is driven by the Lorentz force around the Fermi contour 
in the fc-plane. In the x — y plane, it executes an orbit similar in shape but scaled in dimensions 
by h/\e\B z and turned through tt/2. It has also a z-component of motion which will be discussed 
later. 

The period T of the cyclotron motion, or equivalently the cyclotron frequency lo c , are ob- 
tained from (6,7) by the direct integration. Since the magnitude of these parameters depend lin- 
early on the amplitude of the perpendicular field component, it is more convenient to characterize 
the orbits by the cyclotron effective mass m c defined with a help of the equation to c = \e\B z /m c . 
The explicite expression relating m c to the shape of the Fermi contour reads 

h 2 f dk 

where dk denotes an element of a length of the Fermi line. More often the equivalent expression 
(2) is used. 

It is well known that, unlike in the 3D case, the area Sf surrounded by the /c-space orbit is 
in 2D systems identical to the Fermi area which is determined by the concentration N e of free 
carriers. Since the density of states g is related to the concentration by g = dN e /dE we can 
write 
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where a spin degeneracy is included. Making use of these expressions and of the equation (2) a 
simple relation between the density of states and the cyclotron mass is obtained 



9 = (10) 

The quantization of the in-plane electron motion by B z can be taken into account using the 
Bohr-Sommerfeld quasi-classical quantization rules which state that each quantized real space 
trajectory encloses an integer number of flux quanta /i/|e|. This procedure yields a discrete 
spectrum of Landau levels and the density of states becomes a series of delta functions separated 
by huj c . Note that the filling factor of each level is 2\e\B z /h, i.e. independent of the cyclotron 
mass and the same as in the case of exactly perpendicular magnetic field. 

To estimate the validity of our approximate treatment of the electron spectra we consider a 
model in which a 2D system is confined to the x — y plane by a harmonic potential V con f(z) = 
mil 2 z 2 /2. In this special case the eigenenergies can be found in an analytic form both for 
in-plane and tilted magnetic fields Assuming first, as before, B z = we can write 

1 h 2 k 2 h 2 k 2 

E n (k x ,k y ) = fiuj(n + -) + —^ + —^, n = 0,l,... (11) 

2 2m 2m 

where u = (u 2 + S7 2 ) 1 / 2 , uj y = \e\B y /m and m = mu> 2 /£l 2 . 

The magnetic field induced anisotropy is described by the new effective mass m for the x- 
component of the electron motion which grows up with B y . The Fermi contour is distorted 
from the circular shape to the eliptic one. The zero-field Fermi radius Uqf is determined by 
Eqf = h 2 k 2 F /2m where Eqf is the zero- field Fermi energy. Taking into account the changes 
in the energy spectrum and the Fermi energy due to the in-plane magnetic field, the equation 
Ep = E n (k x , k y ) for determination of the Fermi contour can be rewritten in the form 

^OF = ~ + o (-^) 
uj \l u 

From this equation the cyclotron effective mass can be evaluated using (8) and we get 



m r 



V mm. (13) 



The corresponding cyclotron frequency can be written as u c = Qoj z /u>. 

This approximate results should be compared with the exact solutions obtained for the 
parabolic quantum well subjected to tilted magnetic fields: 



to 



2 + {lo z + Q) 2 T Jul + (w* - 



wi,2 = y 1 • (14) 

It is easy to show that the lowest order expansion of these expressions with respect to oj z yields 
exactly the results obtained by our approximate method 

U Z i w 2 ~ S. (15) 
Both exact and approximate solutions are shown in figure 1. 

3. Self-consistent calculations 

The standard semi-empirical model working quantitatively for the lowest conduction states of 
GaAs / AlGaAs heterostructures is used to solve the Schrodinger equation in the envelope function 
approximation. The envelope function is assumed to be built from host quantum states belonging 
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to a single parabolic band. Since the effect of the effective mass mismatch is completely neglected 
and the envelope functions of GaAs and AlGaAs are smoothly matched at the interface, the 
Schrodinger equation has a form given by (1). 
The confining potential 

V conf (z) = V b (z) + V s . c .(z) (16) 

is a sum of the step function Vb(z) = V&0(— z) corresponding to the conduction band disconti- 
nuity between AlGaAs and GaAs and of a term describing the interaction of an electron with 
ions and the electron-electron interaction. This term should be calculated self-consistently and 
can be written as 

V s . c Xz) = V H {z) + V xc (z). (17) 
The Hartree term Vh is determined from the Poisson equation 

d^VnJM^ (18) 
az l e 

and we use an expression calculated by Ruden and Dohler |l^| in a density-functional formalism 
for the exchange-correlation term V xc : 

e 2 {3N e (z)\ 1/3 , s 

The conduction band offset V& and the dielectric constant e enter our calculations as input 
parameters. 

For modulation doped GaAs/AlGaAs heterostructures the total charge density g{z) in the 
equation (18) can be split into parts corresponding to concentrations of electrons, N e (z), their 
parent donors in AlGaAs, N^(z), and ionized residual acceptors in GaAs, N~(z): 

g(z) = e [N e (z) - N+(z) + N-(z)] . (20) 

We accept a usual approximation of constant impurity concentrations and assume donors and 
acceptors to be ionized within certain finite intervals Id and l a : Nf(z) = for —l^—w < z < —w 
and N~(z) = N a for < z < l a , w is the spacer thickness. 

In our calculation, we consider a GaAs/AlGaAs heterostructure with parameters = 
2 x 10 18 cm~ 3 , N a = 10 14 cm -3 , the band offset V& = 225 meV and the dielectric constant 
e = 12.9. Three selected values of the spacer thickness W\ = 40 nm, W2 = 20 nm and 1^3 = 10 
nm yield electron systems with three different concentrations of electrons ./V e i ~ 1.8 x 10 11 cm -2 , 
N e 2 ~ 3.4 x 10 11 cm -2 and A^ e 3 ~ 5.4 x 10 11 cm -2 , respectively. All these systems have only one 
occupied subband. The parameters l a and Id are determined in the course of the self-consistency 
procedure. For more details see [6]. 

The 'egg-like' Fermi contours calculated for the above three concentration and several mag- 
nitudes of the in-plane magnetic fields are shown in figure 2. The real space trajectories have 
similar shapes. Let us note that, since the average value of the out-of-plane coordinate of an 
electron (z)k x is a function of k x , i.e. the position on the Fermi contour , the resulting trajectory 
does not lie exactly in the x — y plane but is slightly tilted. It is elongated in the y direction 
and an electron is close to the interface at wide end of the trajectory and in the bulk at its nar- 
row end. A schematics of a classical real space trajectory illustrating this behaviour is shown, 
together with its projection to the x — y plane, in figure 3. 

The field dependence of the cyclotron effective mass resulting from the self-consistent cal- 
culations is presented in figure 4. For all three electron concentrations the deviations of the 
cyclotron effective mass reach almost 25% of its zero field value in magnetic field B y = 10 T. 
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4. Discussion and summary 



The electronic structure of a 2D electron gas in GaAs/AlGaAs heterostructure subjected to a 
tilted magnetic field has been studied theoretically. We combine (i) the self-consistent quantum 
mechanical calculations of the electron subbands of the 2D electron system in the presence of the 
parallel magnetic field and (ii) the subsequent quasi-classical description of the in-plane electron 
motion under the influence of the perpendicular magnetic field component. 

The parallel magnetic field combined with the confining effect of the quantum well is the 
reason for deviations of the Fermi contour from a circular shape. Note that in this case the 
nature of changes is different than in case of a standard modification of the Fermi line due to 
the periodic potential. Here, the time-reversal symmetry is broken by the magnetic field and, 
consequently, the Fermi contour has no inversion symmetry in the k space, if the quantum well 
is asymmetric. 

The quasi-classical quantization of the in-plane electron motion yields the eigenenergies linear 
in the perpendicular component of the magnetic field. The method is applicable if the energy 
separation of subbands is larger than the separation of the Landau levels. It does not necessarily 
mean that the parallel magnetic field component must be always greater than the perpendicular 
one, but the in-plane field increases the subband separation. 

The deviations of eigenenergies from the linear dependence for higher perpendicular magnetic 
fields are attributed to the fact that the real space trajectories do not lie exactly in the x — y 
plane, as already mentioned in previous sections. In this case the quantization rules should 
not take into account only the perpendicular component of the magnetic field and the area 
surrounded by the projection of a trajectory to the x — y plane but the full field and the area of 
the trajectory itself. 

The in-plane magnetic field dependencies of the cyclotron mass were calculated in the limit 
of the infinitesimally small perpendicular component of the magnetic field. The changes of the 
cyclotron mass are surprisingly large and, to our knowledge, fully in the scope of the present 
experimental technique. 
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Figure captions 



Figure 1. Eigenfrequencies of a parabolic quantum well subjected to a tilted magnetic field, 
(p = 60°. Dashed lines denote the exact solution, full lines the approximate one. 

Figure 2. Self-consistently calculated Fermi lines corresponding to electron concentrations 
N el « 1.8 x 10 11 cm" 2 , N e2 « 3.4 x 10 11 cm" 2 , N e3 « 5.4 x 10 11 cm" 2 for: (a) B=0 T, (b) B=5 
T, (c) B=10 T. 

Figure 3. Schematics of a classical real space trajectory of an electron in an electric and a 
tilted magnetic field (thick line) and its projection to the x — y plane (thin line). 

Figure 4. Self-consistently calculated in-plane magnetic field dependence of the relative cy- 
clotron mass corresponding to electron concentrations iV e i 1.8 x 10 11 cm -2 , AT e2 rs 3.4 x 10 11 
cm~ 2 , iV e3 « 5.4 x 10 11 cm~ 2 . 



